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Ab initio band theory including correlations due to intra-atomic repulsion is applied to study 
charge disproportionation and charge- and spin-ordering in insulating Nao.sCoC^. Various ordering 
patterns (zigzag and two striped) for four-Co supercells are analyzed before focusing on the observed 
"out-of-phase stripe" pattern of antiferromagnetic Co 4+ spins along charge-ordered stripes. This 
pattern relieves frustration and shows distinct analogies with the cuprate layers: a bipartite lattice of 
antialigned spins, with axes at 90° angles. Substantial distinctions with cuprates are also discussed, 
including the tiny gap of a new variant of "charge transfer" type within the Co 3d system. 

PACS numbers: 71.20.Be,71.27.+a,74.70.-b,75.25.+z 



The Na2;Co02 system, which forms the basis for a 
quasi-two-dimensional transition metal oxide super- 
conductor (T c = 4.5 K) when hydrated,[| shows a 
wide variety of unexplained behaviors in the accessi- 
ble range < x < 1. Countcrintuitivcly (considering 
it is a 2D transition metal oxide) it shows uncorre- 
lated behavior in the normal state for x < 0.50 
(superconductivity arises when x « 0.3 samples are 
hydrated). Then, rather unexpectedly for a hole- 
doped band insulator, it displays correlated behavior 
for x > 0.5 including an enhanced linear specific heat 
coefficient and local moments (Curie- Weiss suscepti- 
bility). Both of these regimes are metallic. Precisely 
at x=0.5, however, it evolves through Na ion order- 
ing, charge ordering, and magnetic ordering transi- 
tions to attain a ground state that is insulating^ 
with a very small gap (few tens of meV). 

Much of the interest in this system lies in the 
triangular arrangement of the Co ions, and the ex- 
pectation that the system should be addressable in 
terms of nonmagnetic Co 3+ and spin-half Co 4+ ions. 
Ordering phenomena, whether charge, spin, or or- 
bital, acquires a different character on a triangular 
lattice |M Q than on the heavily studied square lat- 
tice of the cuprate superconductors. The system be- 
comes magnetically ordered at x > 0.75 (antialigned 
stacking of ferromagnetic layers 0). The insulating 
state was discovered by Foo et ai.Q, who presented 
electron diffraction data indicating robust Na ion or- 
dering (in an orthorhombic four-Co supercellj^, 
whose cell shape is shown in Fig. ^ The Na or- 
dering persisted to above room temperature, and it 
was suggested that Na ordering could be coupled to 
charge (hole) ordering. The zigzag Na order they in- 
ferred, involving equally the two distinct types of Na 
sites, was confirmed by Yang et aL^fji an d calcula- 
tions by Zhang et al. concluded that this ordering 
is favored because it minimizes the Coulomb inter- 
actions between the Na + ions.[Tl| 

As the temperature is lowered, a kink in the in- 



plane susceptibility Xab at T c i=88 K signals antifer- 
romagnetic (AFM) ordering of some Co spins. 0,0, 
fl3j Infrared reflectivity studies0 0, detect a 
gap of ~15 meV opening below T c2 = 52 K, where 
Foo et al. observed 3] the onset of insulating behav- 
ior in the resistivity p(T). T C 2 has been called the 
charge-ordering temperature but there is also addi- 
tional ma gne tic rearrangement, signaled by a kink 
in Xc-[13j liil At T C 3=27 K Gasparovic et al. ob- 
served an additional kink in Xab with no signature 
in p(T); this is the temperature where Foo et al. had 
observed structure in p(T) reflecting more highly in- 
sulating behavior. Unlike the upper two transi- 
tions, there is no entropy change 8] at T c3 . The 
interpretation of this onset of insulating behavior 
was suspected to be charge ordering, pi Till invit- 
ing neutron diffraction studies. The two studies re- 
ported to date confirm that (1) there are two 
types of Co ions (i.e. charge disproportionation), 
one consistent with spin half but with a reduced or- 
dered moment (0.25-0.34 ps), the other spin being 
much smaller, and (2) AFM ordering below T c i=88 
K seems to be the type shown in Fig. ^b). 

First principles local density approximation 
(LDA) calculations by Singh at x=0.5 predict^} 
ferromagnetic (FM) Co layers to be much more sta- 
ble than a simple AFM arrangement, and treat- 
ing the Na ions explicitly does not change this 
conclusion. ^| A crucial development occurred with 
the discovery that, by using the correlated LDA+U 
method, explicit charge disproportionation occursQ 
above a critical value U cr . Subsequent ordering of 
the holes then can relieve the frustration on the tri- 
angular lattice, whereupon it was found that AFM 
order became favored over FM at cc=0.5.[T^ 

The questions of type(s) of order, mechanism of 
ordering, and character of the insulating state have 
begun to be clarified by the data discussed above. A 
quantitative understanding of the behavior is likely 
to require accounting for: the multiband t2 g system 
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FIG. 1: (color online) Charge and spin ordering of 
(a)zigzag (ZZ) and (b) stripe (OP-ST) patterns sug- 
gested by Gasparovic et al. and Yokoi et al., and (c) 
in-phase stripe in the Co layer. In the OP-ST model (b) 
the Co2 ions lies at a site of in-plane inversion symmetry 
and is neighbored symmetrically by | and J. spins. Trian- 
gle and filled circle denote magnetic (Col) and nonmag- 
netic (Co2) cobalts, respectively, and the direction of the 
triangles indicates spin orientation. In the calculations, 
Na lies above Co2. 



with symmetry broken down to singlet a g + doublet 
e' g by the hexagonal ligand field; the triangular, non- 
bipartite lattice that frustrates AFM ordering and 
provides several nearly degenerate possibilities for 
charge order |2(], |2l|; correlation effects strong 
enough to drive charge disproportionation but small 
enough to leave a tiny charge gap. 

We have addressed these questions using the same 
methods as for our previous studies on this 
system, |2j, [19| with attention to Brillouin zone sam- 
pling (up to 312 k-points in the irreducible zone). 
Specifically, we utilize the disproportionated states 
provided by the LDA+U approach to address, for 
the three superlattice symmetries shown in Fig. ^ 
the energetics, the relative orientations and magni- 
tudes of the magnetic moments, and characteristics 
of the electronic structure in the insulating phase. 
These Co orderings correspond to (a) zigzag (ZZ), 
(b) out-of-phase stripe (OP-ST), and (c) in-phase 
stripe (IP-ST). Note that the Na ion zigzag phase is 
not the same as this Co ZZ order; the Na zigzag in- 
volves one site on top of Co and another site not on 
top of any Co that is "less zigzag" than this Co ZZ 
order. There is some controversy of the effect of Na 
order: Na-Na repulsion accounts for the preferred or- 
der |ll| without further mechanisms; Na + order and 
Co 3+ /Co 4+ order are coupled [23J; or the observed 
Co 3+ /Co 4+ order can be obtained from a single- 
band extended Hubbard model 12 ll without any fur- 
ther complications. Li et al. 18] used the observed 
Na superstructure |8( and included O ion relaxation 
but did not compare the results with simpler Na ar- 
rangements [24|. We do not address this question 
here but note that the Na ordering we adopt, used 
in previous work,^ij serves to break the the sym- 



0.7 
0.6 
0.5 



~«0.4 
0.3 



0.1 




JET 



XT 



.1 
0.2 



B" El St Col 

0—0 st Co2 

A A ZZ Col 

v- — v zz Co2 



X — 6 06 — * — ' — — g — s 



2 3 
U(eV) 



FIG. 2: (color online)Effect of U on magnitude of the 
Co local magnetic moments rm in the stripe and zigzag 
patterns of AFM Nao.sCo0 2 . At U c =1.5 eV, charge 
disproportionation Col^Co 4+ (S=l/2) and Co2^Co 3+ 
(S=0), occurs with gap opening. 



metry of the Co sites. Specifically, the Na ions sit 
above the Co2 site (which becomes the nonmagnetic 
Co 3 + site). 

Our attempts to obtain an AFM state for the IP- 
ST model of Fig. [Jc) converged to a FM state (or 
nonmagnetic, if AFM symmetry was enforced). |25j| 
Thus we will consider only the model given in Fig. 
njb) as a stripe (ST) pattern in this paper. Within 
LDA, FM is favored over AFM for both ZZ and ST 
patterns, as for all other values of x. The energy 
difference is substantial for ZZ (300 meV/Co) but 
surprisingly small for ST (8 meV/Co). This favoring 
of FM by LDA confirms the need for including effects 
of correlation, as we do by applying the LDA+U 
method. As emphasized previously, results depend 
on the value of U and it is necessary to determine the 
appropriate value. This is rather straightforward for 
£=0.5: there should be disproportionation, charge- 
and spin-ordering (AFM) and a very small gap. We 
first review behavior versus the repulsion strength 
U. Hund's rule J=l eV is kept fixed. 

The effect of U is evident in the calculated Co mo- 
ments, displayed in Fig. |2 At U=0 (i.e. LDA level) 
the effects of symmetry (determined by the Na place- 
ment) is already strong. The Col and Co2 moments 
are nearly equal for ZZ, while there is already almost 
negligible moment on Co2 for ST. This difference re- 
flects the higher symmetry of the Co2 ion in the ST 
pattern of Fig. ^b) : it is surrounded symmetrically 
(in-plane inversion) by two Coif, two ColJ., and two 
nonmagnetic Co2. Increasing U, Col magnetic mo- 
ments in both patterns increase monotonically and 
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become identical at and above £7=1.5 eV. The low- 
spin Co2 magnetic moments show a much greater 
difference between the two patterns. For ST it be- 
comes immediately (by £7=0.5 eV) zero, while for 
ZZ there is at £7=0.5 eV what might be identified as 
the charge disproportionation transition, but beyond 
this point the Co2 moment simply decreases mono- 
tonically, never becoming identically zero. From Fig. 
^a) the lack of symmetry in ZZ is clear: although 
surrounded by two | and two J, Col spins, and two 
low-spin Co2 ions, there is no in-plane inversion, so 
a moment is allowed. While the Col and Co2 ions 
are clearly disproportionated in Fig. [21 the charge 
difference is only ^0.2 e (this difference is 0.02 e 
smaller for ST than ZZ). 

The gap opens, for both patterns, at £7=£7 cr =1.5 
eV. It is noteworthy that the disproportionation had 
occurred already at smaller U (see Fig. [2J; thus 
we find here a richer behavior than in our previous 
studies with smaller cells or different Na concentra- 
tions, where disproportionation/ordering had coin- 
cided with gap opening. Such a difference was also 
obtained for a similar supercell by Li et al, who 
however found a somewhat larger value of U^.[l8| 
This critical interaction strength coincides with a 
Mott-like transition in the Col a g states, with up- 
per and lower Hubbard bands separated by 2.2 eV 
as shown in Fig. |2| The distinctive a g state is 
prominent already within LDA, arising from sym- 
metry breaking due to the hexagonal ligand field 
[t2 g — > a g +e p g rime), and the gap opening at U cr con- 
stitutes the band description of an orbitally selective 
Mott transition. Almost independent of £7, the ZZ 
pattern is favored over ST by the very small value 
of 22 meV/Co; given the fact that the Na zigzag ar- 
rangement has been simplified in our calculations, 
we can conclude that these two patterns have no 
significant difference in energy. 

The upper Hubbard band in the insulating state 
is quite flat, and in fact shows more dispersion per- 
pendicular to the layers (200 meV) than within the 
layers. The band structures are pictured along di- 
rections parallel to (x), and perpendicular to (y), the 
ZZ or ST chains of Co ions in Fig. for £7=2 eV, 
a value slightly above U cr to make the gaps more 
clearly visible. For the ZZ case, the minimum gap 
occurs at a corner of the zone that is not shown. 
The dispersions going away from the zone center of 
the uppermost two valence bands are entirely differ- 
ent for the two patterns, being positive for ST but 
negative for ZZ. The uppermost bands have primar- 
ily Co2 e' g (t2 g ), and not a g , character. In spite of 
these differences in dispersion through the zone, the 
orbital-projected density of states, not shown, is ex- 
tremely similar for the ZZ and ST bands shown in 
Fig. El 
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FIG. 3: (color online) View of the AFM band structures 
in the ti g manifold at U=2 eV for (a) the zigzag (ZZ) and 
(b)the stripe (ST) patterns. The plot is along perpen- 
dicular (r-Y) and parallel (r-X) directions for the each 
chain. Z denotes the zone-boundary point along < 001 > 
direction. The thickened lines emphasize the bands with 
strong a g character for each spin of a magnetic Col. The 
energy zero lies in the gap. 

We focus now on the observed ST pattern. The 
gap occurs at the zone boundary point Y along k y . 
Although Co 3d states border the gap on both sides, 
this is an unusual d — d charge transfer gap (not 
the usual p — d case), with unoccupied Col minority 
a g states above, and primarily Co2 e' g states below 
the gap. The stronger dispersion of the upper Hub- 
bard band above the gap along Y-r compared to X~r 
can be understood as follows. Electrons excited into 
the upper Hubbard reside in the minority a g states 
on Col, for example, a spin f electron will hop be- 
tween Col ions with moments oriented J.. Propagat- 
ing in the y direction, it can hop through a single 
Co2 ion; in the x direction however, it must avoid 
the £7 cost of hopping onto an oppositely aligned 
Col ion, thus requiring hops through two Co2 ions 
before returning to another Col J, ion, and its dis- 
persion is reduced accordingly. Valence band holes 
introduced into the system will occupy nonmagnetic 
Co2 e' g states, while electrons will occupy minority 
a g states on Col. 

Since the crystal field (t 2g — > e g ) gap is ~2 eV, 
the optical transitions in the IR for the magnetically 
disordered metallic phase (T>T c i) reflect e' g — * a g 
excitations, i.e. transitions within the t 2g complex. 
Below the metal- insulator transition at T C 2, the ex- 
citations across the gap are to the upper Hubbard 
band, and the main weight of these transitions - the 
new (Col «-> Co2) charge-transfer type mentioned 
above - is shifted up in energy by only a few tens of 
meV. EE El This small shift is consistent with 
the small bandwidth that we find for the upper Hub- 
bard band. 

An interesting analogy with the AFM cuprate 
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layer arises. In this observed ST pattern Fig. ^b), 
the antiferromagnetic arrangement has the bipar- 
tite, and 90°, topology of spins characteristic of 
the cuprate plane, except with anisotropy of (Co + - 
Co 4+ ) parallel and perpendicular hopping ampli- 
tudes t x ,t y and exchange couplings J x , J y . Note that 
exchange J x is between near neighbors while J y is 
between second Co neighbors. There are however 
strong distinctions to be made with cuprates. In 
cuprates U is 3-4 times larger, the effective metal- 
metal near neighbor hopping is 2-3 times greater, 
and the 2 eV gap typifies a robust Mott insula- 
tor. In this cobaltate system the tiny gap reflects 
a marginally insulating correlated state, and the 
low energy excitations require three d bands, ver- 
sus the dominance of the single d x i_ y i state in the 
cuprates. It also seems that magnetic coupling can- 
not be treated in the usual Heisenberg form, because 
the superexchange mechanism is not dominant and 
the rather flimsy moments depend strongly on the 
type of magnetic order. (A spin Hamiltonian might 
be reasonable for treatment of spin waves within a 
given ordered state.) 

Our results suggest a specific picture of the tem- 
perature evolution at x=0.5. Noting that the FM 
ordered layers for x >0.75 are consistent with itin- 
erant character and the 0.5 < x < 0.75 regime with 
fluctuation-suppressed magnetism, the magnetic or- 
dering below T c i may be more of a spin density wave 
(SDW) character which gaps some but not all of the 
Fermi surface; recall that for the ST pattern the FM- 
AFM energy difference is very small at small U, and 
that even at U=0 there is a substantial difference in 
moments on Col and Co2, that is, an SDW . Several 
band structure studies have pointed out nesting fea- 
tures in the paramagnetic Fermi surface. 0, 0, |2(| 
The challenge that this picture must face is that the 
primary magnetic order is unchanged at the insu- 
lating transition T C 2~. Gasparovic et al. find that 



the principal ordered moment grows with decreas- 
ing temperature 13] continuously through the insu- 
lating transition at T C 2. The additional order that 
results in a kink in \c at T C 2 has not yet been elu- 
cidated, but our results are consistent with the pre- 
vailing picture that disproportionation arises finally 
at T C 2. The redistribution of spectral weight below 
T C 2 observed in optical experiments show differences 
(weight shifted to 20-30 meVpjl or 70-100 meV0), 
but they seem consistent with the correlated band 
structure of Fig. [3] and particularly the narrowness 
of the unoccupied band. 

It can reasonably be asked whether the ground 
state of this system should be considered as a cor- 
related insulator, as outlined above, or instead as 
perhaps a SDW (at T ci ) - CDW (at T c2 ) system. 
Balicas et al. have found that, when an applied 
in-plane field increases beyond 25 T the conductiv- 
ity increases by a factor of two (a sort of insulator- 
metal transition), and observation of magnetoresis- 
tance oscillations sug gest the restoration of part of 
the Fermi surface. |27j Certainly the insulating phase 
is delicate. However, the observation of a substan- 
tial ordered moment on Col (0.25-0.34 ur)[HII^| 
and little or none on Co2 speaks for a charge dis- 
proportion picture (into identifiable Co 4+ and Co 3+ 
moments) and hence a correlated insulator below 
T C 2. The observed value of the Col moment is re- 
duced somewhat from our calculated value of 0.5/is, 
as would be expected from two-dimensional fluctu- 
ations of a spin-half moment already reduced sub- 
stantially by hvbridizationp^ with O 2p states. 
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